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ABSTRACT 


An  expansion  of  a  formally  rigorous  integral  solution  of  the 

-  / ,  .  ,-y 

Boltzmann  equation  in  powers  of  the  electric  field  and  used  to 

demonstrate  clearly  the  correct  expression  for  the  conductivity  of  a 
spatially  homogeneous  slightly  ionized  gas.  The  conditions  under  which 
the  isotropic  part  of  the  electron  distribution  function  reduces  to  the 
Margenau  distribution  are  determined.  A  discussion  of  the  physical 
reasons  for  the  possibility  of  negative  conductivity  is  given.  . 


I.  INTRODUCTION 


The  theories  of  the  interaction  of  slightly  ionized  gases 

with  electromagnetic  fields  as  developed  by  many  authors  2*  3>  4’  5* 

7 

are  usually  based  on  a  modified  Lorentz  gas  approximation.  That  is, 
electron-electron  collisions  as  well  as  ion-ion  collisions  are  ignored, 

.‘Vi 

and  the  recoil  of  the  neutral  molecules  is  taken  into  account  to  first 
order  only.  In  addition  the  Boltzmann  equation  is  expanded  in  Legendre 
polynomials  in  velocity  space.  The  latter  procedure  leads  to  an  infinite 
set  nf  coupled  differential  equations  for  the  expansion  coefficients  which 
is  truncated,  usually  by  setting  all  but  the  first  two  coefficients  equal 
to  zero.  This  procedure  leads  to  the  following  expression  for  the  cur¬ 
rent  density: 


4tre  /»*  .  ‘3  *f0(v)dv 

3  J0  X>  *v 


(1) 


in  which  ( - e)  is  the  charge  of  an  electron  and  f0(v)  is  the  isotropic  part 
of  the  electron  distribution  function.  For  an  electric  field  of  the  form 


<(t)  =  Ecoseit 


(2) 


1 


the  quantity  <  x  >  is  given  by  the  equation 


•  -eh, 

<  x>  =  - - - —  vcoswt  +  wsinuit  1  (3) 

m(vz  +  tan  l-  J 

where  v  is  the  electron  collision  frequency  (normally  a  function  of 
velocity). 

g 

In  a  later  paper,  Margenau  has  shown  that  equation  (1)  is 
valid  for  more  general  conditions  than  initially  had  been  demonstrated. 

Q 

It  was  then  pointed  out  that  equation  (1)  does  not  in  itself  guarantee  the 
positive  definiteness  of  the  electric  conductivity,  and  that,  in  fact,  sub¬ 
stitution  of  a  8 -function  in  the  magnitude  of  the  electron  velocity  for  fQ(v) 
and  the  experimental  collision  frequency  in  air  for  v  ,  yields  a  negative 
value  for  the  conductivity.  The  following,  alternate  formula  for  the  current 

a 

density  was  proposed: 

j  =  -4*6^'  <x>  v2f0(v)dv.  (4) 

In  order  to  obtain  this  formula,  which  guarantees  positive 
conductivities,  the  philosophy  was  adopted  that  the  only  effect  of  the  electric 
field  (which  we  take  in  the  x  direction  for  definiteness)  is  to  shift  the 
"equilibrium  distribution",  fQ(v)  in  the  following  fashion: 


2 


(5) 


f(v 


x’  y* 


t) 


fo(vx- 


y  v2 


)  (l  .  _f£_)  , 

*Vv 


where  f(v,  t)  is  the  new  distribution  function,  g  is  the  velocity  shift  of 

an  average  electron,  and  the  factor  (1  -  — — )  accounts  for  the  fact 

»vx 

that  the  volume  element  in  the  "shifted"  region  differs  from  the  unshifted 
volume  element,  g  was  shown  to  be  equal  to  <  x>  as  given  by  equation  (3). 

When  the  current  density  is  calculated  using  a  distribution 

» 

function  of  the  form  given  by  equation  (5)  and  is  subsequently  expanded 
in  powers  of  <  k  >  ,  the  lowest  order  non-vanishing  result  is  given  by 

equation  (4). 

Taylor1  °has  recently  computed  the  average  shift  of  the  electron 
velocity  bv  what  is  essentially  a  random  walk  technique.  His  result  agrees 
with  equation  (3).  Subsequently  he  uses  equation  (5)  for  the  calculation 
of  the  current  density  to  all  orders  in  <  x>  and,  as  a  special  case,  ob¬ 
tains  equation  (4)  to  first  order  in  <x>  ,  Taylor's  calculation  of  <x>  , 

however,  is  correct  only  if  the  collision  frequency  is  independent  of  velo¬ 
city,  because  he  neglects  to  take  into  account  certain  retardation  effects, 
which  we  shall  have  occasion  to  discuss  later  on  in  this  paper  (equation  ( 15)  ). 

In  Section  II  of  this  paper  we  expand  a  formally  rigorous  solu¬ 
tion  of  the  Boltzmann  equation  for  a  modified  Lorentz  gas  in  powers  of 
the  electric  field  and  in  powers  of  the  ratio  of  the  electron  mass  to  the 
mass  of  the  molecule.  The  function  f0(v,  t)  is  precisely  defined  and  the 
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r 


conditions  under  which  it  reduces  to  the  steady  state  Margenau  distri¬ 
bution^  are  determined.  In  Section  III  the  current  density  is  computed 
and  shown  to  be  given  by  equation  ( 1)  to  first  order  in  the  electric  field. 
The  physical  reasons  for  the  possibility  of  negative  conductivities  are 
discussed. 


'<*W 


II.  THE  ELECTRON  DISTRIBUTION  FUNCTION 

The  Boltzmann  equation  for  a  modified  Lorentz  gas  which 
is  homogeneous  in  space  and  on  which  is  imposed  an  electric  field  of  the 
form  given  by  equation  (2)  is 


*f(v.  t) 

•  t 


+  ycoswt 


#f{v,  t) 


=  -  v(v)f(v,  t)  +  S(y,  t) 


where  Y=  — "  ■  E  and  S(v,  t)  is  a  source  term.  The  collision  frequency 
v (v)  describes  all  types  of  collisions  -  clastic,  inelastic,  ionizing  colli¬ 
sions,  absorption,  etc.  S(_v,  t)  is  the  rate  at  which  the  number  of  electrons 
in  the  velocity  range  v,  v  4  dy  is  added  due  to  elastic  collisions  as  well 
as  due  to  all  other  processes.  In  general  S(v,  t)  depends  on  the  form  of 


f(v,  t). 


To  obtain  a  formal  solution  of  equation  (6)  we  first  make 


the  following  transformation: 


sinwt  =  c  , 


Let  F(c,  Vy,  vz,  t)  =  f(vx,  Vy,  vz,  t),  then  the  equation  becomes 


—  3  in  (at,  v, 
0)  ’ 

(S) 


V*  ''t*  ^  y 

•  +  v(c  +  -^j-sinut,  vy,  vz)F(c,  vy,  vz>  t)  =  S(c  +_ 


The  solution  of  this  equation  is: 


t)  =  TsCe +-^~  sin  t’ .  vy,  vz,  t')exp£-^* v(c  +  sin  t",  vyl  vz)dt"]d1 
*o  t* 


+  F(c,  vy,  Vz,  t0)exp  v(c  +  -jj-  sin  t\  vy,  vz)dt'j  , 


(9) 


or,  after  transformation  to  the  original  variables: 


t  t 

t)-J  S(vx  -  x(t,  t’),  vy,  Vz,  t1)  exp  [  -  J  V(vx  -  -^-x(t,  t”),  vyl  vz)dt"]dt' 


(10) 


+  exp  [-/  V(vx  -  x(t,  t'),  vy,  vz,  t')dt'j  f(vx  -  x(t,  tD),  vy,  Vz,  tD) 


where  the  notation 


x(t,  t')  =  [sinwt  -  sinwt]  (11) 

has  been  introduced,  and  where  t0  is  some  initial  time  at  which  f  is 
specified.  Equation  (10)  is,  of  course,  not  a  solution  of  the  problem  since 
S(v,  t)  is  itself  related  to  f(v,  t).  It  is  merely  an  integral  form  of  equation 
(6)  and  is  completely  equivalent  to  it.  Equation  (10)  does,  however  provide 
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a  convenient  starting  point  for  a  systematic  computation  of  the  distribu¬ 
tion  function  as  a  series  in  powers  of  Y . 

Consider  the  case  where  t  s  -  «  .  Because  of  the  posi¬ 

tive  definiteness  of  the  collision  frequency,  v,  the  second  term  on  the 
right  had  side  vanishes  for  all  t  >  -  oo  ,  leading  to  the  equation 

t  t 

f<vx,  v  .  v  t) -f  S(vx  -  x(t,  t'),  v  .  vz,  t1)  exp  v  T  vx  -  x(t.  t"),  vy>  vzjdt"ldt'  . 

'  1  *  -00  t' 

(12) 

We  digress  temporarily  to  discuss  certain  aspects  of  Taylor's 
paper.  Assume  that  the  molecules  are  infinitely  heavy  and  that  the  scat¬ 
tering  cross  section  is  isotropic.  In  this  case 


S(v,  t)  =  N 


v'f(v|,  t)  «r ( vj 


(13) 


=  v(v)fQ(v,  t)  =  S0(v,  t) 


where  fQ( v,  t)  and  S0(v,  t)  are  defined  by  the  equations 

fo(v*  t]  =  *>dn  5  s0(v'  =  s<v’  t,d0  •  (14> 

Here  N  is  the  number  of  molecules  per  unit  volume  and  »(v.  -svldfl 
is  the  differential  scattering  cross  section.  The  final  form  of  equation  (13) 
results  because  in  this  particular  case  v1  =  v  and  0  is  a  function  of  the 
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? 


velocity  magnitude  alone.  We  now  multiply  equation  (12)  by  v  ,  integrate 

•  X 

over  all  velocities  and  divide  by  the  electron  density  n  to  obtain  the  average 
electron  velocity.  After  substitution  for  S(v,  t)  from  equation  (13),  we 
obtain 


<  V  =  4  If  \v{\  -  *(t-  v)'  v  v  v)  fo(vx 


-  eo 


x(t,  t’),  V  , 

y 


-  X(t,  t"),  vy, 


dt’  d 


3 

V  . 


v  ,  f) 
z 

(15) 


For  the  special  case  of  velocity  independent  collision  frequency  and  time 
independent  fQ.  equation  (15)  reduces  to 


<  v 


X 


x(t, 


t ')  V  e 


-  V  (t  -  V) 


dt' 


(16) 


which  upon  slight  change  in  notation  reduces  to  equation  (5)  of  Taylor's 
paper,  and  the  current  density  is  a  linear  function  of  the  electric  field 
regardless  of  the  strength  of  the  field.  In  general,  as  can  be  seen,  a 
velocity  dependence  of  v  implies  a  complicated  time  dependence  of  V 
in  integrals  such  as  (16).  If  equation  (15)  is  expanded  in  powers  of  y 
(again,  with  isotropic,  time  independent  fQ)  <vx>  1*  not  given  by  equa¬ 
tion  (16)  but  contains  an  additional  term  of  the  same  order,  which  is  pro¬ 
portional  to  — — -  ,  and  higher  order  terms  in  Y.  This  equation  is 

»vx 

therefore  not  adequate  for  the  calculation  of  the  conductivity  even  to  lowest 
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order  non- vanishing  terms  in  Y,  much  less  for  the  computation  of 
higher  order  terms  as  Taylor  has  done. 

We  now  return  to  the  computation  of  the  distribution  function. 
We  expand  equation  (12)  in  powers  of  Y  to  obtain 

f(v,  t)  -J  e'V(t  '  ‘^dt'^v,  f)  -  -^-[a^t,  t')  *  a2(t*  t')V"  1  S(v,  t’)] 

Y2  f r  a2<‘*  *’>  av  1  2 

+  ^tU— —7b\  s<* 


2b)  ‘ 


f) 


t')  y 


(17> 


2al(t’  t,,a2(t’  0  8  V  8 S(v,  t')  ,  .2 . #2S(v,  t')  )\ 

5  st  »vx  +ai(t,t)T^2 — r;  . 


where  V  =  v(v)  and 


a,(t,  t')  =  i(t,  t')  ;a2  *(t,  t")dt' 


a32(t,  t')  =  f  x2{t,  t 

Y  t' 


(18) 


We  now  integrate  both  sides  of  equation  (17)  over  solid  angle  in  velocity 
space,  and  obtain,  under  the  assumption  that  S(v,  t)  is  isotropic,  the 
following  equation: 


With  the  exception  of  a  slight  digression  to  a  discussion  of  more  general 
conditions  following  equation  (49)  in  Section  111,  we  shall  restrict  our¬ 
selves  to  the  case  where  elastic  scattering  of  electrons  by  neutral  mole¬ 
cules  is  the  only  contribution  to  the  source  term  for  the  remainder  of  this 
paper.  Also,  with  the  exception  of  the  appendix,  we  consider  only  the 
case  where  the  differential  scattering  cross  section  is  isotropic  in  the 
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center  of  mass  system.  Note  that  for  infinitely  heavy  molecules,  these 
assumptions  imply  that  S(v,  t)  has  the  isotropic  value  given  by  equa¬ 
tion  (13).  Substitution  from  this  relationship  into  equation  (20)  shows 

J  fD  2 

that  —  —  is  small  of  order  V  . 

We  now  examine  the  corrections  to  equation  (13)  which 
arise  when  recoil  is  taken  into  account  to  first  order. 

Let  <r(v',  $)  d fi  be  the  differential  cross  section  in  the 
laboratory  system  for  scattering  electrons  through  an  angle  $  into  the 
velocity  range  v,  v  +  d^v.  We  assume  that  the  temperature  of  the  neutral 
molecules  is  low  enough  so  that  their  velocity  before  the  collision  can 
be  ignored;  thus  v1  is  the  velocity  of  the  electron  before  the  collision 
and  0  is  the  scattering  angle  in  the  laboratory  frame  of  reference.  As 
a  consequence  of  the  invariance  of  the  differential  scattering  cross 
section  to  a  transformation  to  the  center  of  mass  system  we  have 

<MV,  0 ) an  =  ac(g,  0c)dnc  (21) 

in  which  the  subscript  "c"  refers  to  center  of  mass  system  and  g  is 
the  magnitude  of  the  relative  velocity.  Using  the  relations  resulting 
from  conservation  of  energy  and  momentum  and  the  fact  that  the  velocity 
of  the  molecule  before  collision  is  taken  as  zero,  we  get  the  equation 


n 


where  m  and  M  are  the  electron  and  molecule  masses,  respectively. 
The  azimuthal  angle  is  invariant  to  the  transformation.  Thus 


d(cos  6  ) 

a(v',  6 )  =  a(g,  $) - -  -  (1  +  2—  cos0)  a  (v\  6),  (23) 

c  c  d(cos  0)  M  c  c 

where  we  used  the  relation  v*  =  g  in  the  last  expression.  In  general, 
assuming  ac(g,  &c)  is  known,  one  would  also  use  equation  (22)  to  substi¬ 
tute  for  0c<  In  the  case  considered  here  the  cross  section  is  isotropic 
in  the  center  of  mass  system.  Thus 

o(v\  0)  =  (1  +  2  —  COS0)  cr(v’)  +  0(-2L)2  .  (24) 

M  c  M 

Following  Morse,  Allis,  and  Lamar2,  we  get  the  following  form  for  the 
source  function: 

S(v,  t)  =  -JL-  f  f(v',  t)  (v')4  ff(v',  0)dfl  (25) 

v3  ‘'n  ~ 

where  N  iB  the  number  of  molecules  per  unit  volume.  In  the  integration 
of  equation  (25)  over  solid  angle  ;v  rather  than  v1  is  held  constant.  Thus 
one  must  express  v1  in  terms  of  v  and  6.  The  required  relation  is 


easily  obtained  by  using  the  equations  expressing  conservation  of  momen¬ 
tum  and  energy.  The  result  is 


v' 


m 

M 


( 1  -  cos 


+  Q(-S— )2  . 

M 


(26) 


We  now  split  the  distribution  function  into  isotropic  and  non-isotropic 
parts, 


f(v,  t)  =  fo(v,  t)  4  fj(v,  t) 


(27) 


and  initially  make  no  assumptions  as  to  the  relative  size  of  the  two  parts. 
Using  equations  (24),  (26)  and  (27)  in  equation  (25)  and  expanding  the  re¬ 
sult  in  powers  of  m/M  we  obtain 


S(v,  t)  =  v(v)f  (v,  t)  4  -2— 

Mv5 


+  — ~  V(v)f  l(v,  t)  -  -2L 
M  1  "  Mv 


i-  £  v3v(v)fo(v,  t)J 

-2-  r  v3  V(v)f  '(v.  t)1  4  0(-2L)2 
8  v  L  1  ”  J  M 


(28) 


where 


■f  11<V,  t)  =  /*  ((v1,  tlcoBjdfl,  v(v)  =N  f  v<r(v,  0  )dfl  =  4irNvff,(v)  +  0(_2-)2  . 

*/n  ~  ■'n  ™ 


(29) 
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Consider  the  integral  expression  for  f  j'  . 


It  can  be  written  as 


and,  since  v'  =  v  +  0(_2L)  one  can  write 


f,1  =  f  f( v * ,  t ) v *  •  v  .-ilL'-l-il.  dV  +  0(J2L)  . 

1  J  ~  ~  ~  v4  M 

Clearly  f  ^(v,  t)  is  an  odd  function  of  the  components  of  v  to  order  -12-  ,  and  is 
1  —  ~  M 

hence,  in  this  approximation,  purely  anisotropic,  that  is 


(30) 


When  S(jv,  t),  as  given  by  equation  (28)  is  substituted  on  the  right  hand 
side  of  equation  (17)  it  is  seen  that  f(v,  t)  is  isotropic  if  terms  of  order 
Y  and  higher  order  terms  arc  negtoctod.  That  is  f](jv,  t)  is  of  order  Y. 
By  making  use  of  this  fact  we  can  rewrite  equation  (28)  as  follows. 


14 


We  can  thus  see  that  when  the  center  of  mass  differential  scattering  cross 

section  is  isotropic,  the  source  term  is  isotropic  through  order  --n  , 

M 

In  the  derivation  of  equation  (19)  and  (20)  we  have  assumed  that 

the  source  is  isotropic.  By  substituting  S(v,  t)  as  given  by  equation  (31) 

into  equation  (17)  and  repeating  the  ensuing  calculation,  (without  the  explicit 

assumption  that  S(v,  t)  is  purely  isotropic),  we  see  that  equations  (19)  and 

2 

(20)  are  still  satisfied  if  terms  of  order  Y  -  and  higher  order  terms 

M 

are  neglected.  In  particular  equation  (20)  takes  on  the  following  form: 


af0(v,  t) 

at 


Mv 


[^qK-]  +  — ^~T  f 

av  u  o  j  6u>  *“ 


e“  V  (t  -  t')  v  dt' 


'  (Vv-  *■>  {  [twY  .)] 

*  7" ( r  17  *  )  i r  [ **“•  *•> *''] } 

*  T  17  17^’  £ 

+  r/i.  _2_+  -®L\f  (v,  t')]  (  °&1  (tl  1 1  )  +o(  y2-2L)  +  o(  y4) 
L\ v  av  #v2 /  °  J  \at  /  m 


(32) 


Equation  (32)  implies  that  f0(v,  t)  is  a  rather  weak  function  of 

time  (  - *°  is  of  order  JZL  or  y2,  whichever  is  larger).  This  fact, 
at  M  1 

together  with  the  fact  that  the  presence  of  the  exponential  e“  v(^"t')  8Up. 
presses  contributions  with  t'<<  t  in  the  integral  dn  the  right  hand  side 
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term  proportional  to  Y  £  we  get 


gfp(v.  t) 

•  t 


Mv 


dv  [  v3v  fo(v*  l)] 


3v 


9  r  2  cj  . 

— “  (cos  ci>  t  +  — - 7-  sinoit  cosci) t) 

9  v  L  v 


2  2 

W  +  V 


>fQ(v>  t) 

9  v 


]  +  0(  y 


(35) 

2  m 
M 


)  +  0(  y4) 


By  making  use  of  the  weak  time  dependence  of  fQ(v,  t),  we  can  integrate 
over  a  period  of  the  electric  field  holding  fQ(v,  t)  constant.  Furthermore, 

if  we  add  to  the  term  proportional  to  m/M  the  well-known  correction  due 

4  5  6 

to  the  thermal  motion  of  the  molecules,  ’  *  we  get  the  steady  state  result 


Y2 

d 

V2  V  dfo(v)  1 

-  m  d 

J  v^V  f  (v)  +  kT 

df0(vn 

6 

dv 

2  2  dv 

L  <0  - f  v  QV  J 

M  dv 

1  v  v  + - 

1  L  mv 

dv  J 

(36) 


which  upon  integration  leads  to  the  Margenau  distribution  function 


fo<v>  =  f0<v!  =  Cexp  f  -/  - mJdV - - - —  1  .  <37> 

*■  o  kT  +  y2M/6(  w2  +  V2)  •* 

in  which  C  is  a  normalization  constant. 

Once  again  we  make  use  of  the  result  that  the  time  dependence  of 
f0(v,  t)  is  insignificant  and  integrate  equation  (33)  over  time,  holding  fQ(v,  t) 
fixed  on  the  right  hand  side.  The  result  can  be  used  to  express  fQ(v,  t') 
in  terms  of  f0(v,  t). 
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(.[v,  t')  =  IJv,  t)  *  (f  -  t)  -2 - 1~  r  V3VUV,  til 

Mv2  »”  L  1 

•t4  llfVH  (TT7t5)f“K,‘ 
*[/W](^^)} .  (3S) 

Now,  by  using  this  equation  together  with  equation  (31)  we  obtain 

s,v,  «■)  ■  vyv.  .>♦[.-  v,,-..,]  -a.  _ L.  [,»vy„  .)]  --££{  }. 

(39) 

where  the  curly  brackets  are  meant  to  contain  the  same  quantity  as  they 
do  in  equation  (38),  We  substitute  the  r.  h.  s.  of  equation  (39)  for  S(v,  t1) 
in  the  zeroth  order  term  of  equation  (17),  while  the  value  Vf0(v,  t)  is  used 
for  S(v,  t1)  in  the  higher  order  terms.  The  term  proportional  to  m/M  and 
the  terms  proportional  to  9v/9vx,  JV/»vx)2,  or  (  9  2  V  /  9  vx2)  times 

f0(v,  t)  drop  out  in  the  integration  over  t  and  we  get 


«V.  .1 .  y  v.  .1  -  JL  b«)  -li-  -iMal 

^  °  CJ  v  Qv 


s  (^T-  -  -f)  { ^  (•  4-  -77  *  -£) 

HO) 

2B2(t)  dy  9f0(v,  t)  x  m  3 

-T--3TT? - )  401  Y-S-' 
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o 


where 


b(t)  -  f  e*  ^  v  dt'a.(t,  t')  = - — — —  ("1/0080)1+  O)sino>tl 

1  v2  +  o.2  L  J 

Bj(t)  =/*  e-V^-t'Jtfdf  [  a2(t,  t')l  =  — — -  [1  +  3a)V8in2o)t  +  (v2.2o,2)co_820)t1 

-»  L  -1  <v2  +  «D2)  L  ( v2  +  4 a) 2)  J* 


and 


(41) 


S2(t)  =/  e*V(t*t')  Vdf  ^a2(t,  f)  -a^t,  t')a2(t,  t')^j 

B  til2  J  O)2  -  V2  (2  ti)3v  -  4  v3o>)8in2o>t  +  (5q>2»2-  V4)cos2o)t  | 
2(l/2  +  W2)  L  O)2  +  V2  («2  +  v 2)  ( v2  +  4o>2)  ^ 


In  obtaining  equation  (40)  we  made  use  of  the  fact  that  Bj(t)  and  B2<t)  can 
be  expressed  as  follows; 

Bl<t)  •/  e'',(t_t’)v2dtl  [y  A j(t")dt"]  ,  B  (t)  =  J  e“V(t't',v^f[/  A  (t")dt"  ]  . 

(42) 

We  see  that  our  solution  for  the  electron  distribution  written  as  an  expansion 
in  powers  of  the  electric  field  turns  out  also  to  be  an  expansion  in  Legendre 
polynomials  in  the  cosine  of  the  angle  between  v  and  the  electric  fiel^  to 
the  order  for  which  we  have  determined  the  coefficients  of  the  expansion. 

In  the  conventional  treatment  the  distribution  function  is  expanded  in  spherical 
harmonics  at  the  outset  and  the  assumption  is  made  that  the  contributions 
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proportional  to  Legendre  polynomials  of  second  order  and  higher  are 
negligible.  The  justification  for  this  assumption  being  that  when  it  is  made 
the  series  appears  to  converge  rapidly  on  the  basis  of  the  two  or  three 
terms  which  are  computed. 

Although  our  method  is  lengthy,  it  seems  to  us  to  be  more  satis¬ 
factory  because  we  have  shown  that  the  contributions  not  explicitly  given  in 
equation  (40)  are  indeed  of  higher  order  in  terms  of  certain  expansion  para¬ 
meters  (  y  and  m/M)  than  those  retained.  Furthermore,  we  have  explicitly 
demonstrated  that  the  oscillating  portion  is  indeed  small. 
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III.  ELECTRICAL  CONDUCTIVITY 


If  we  substitute  the  result  given  by  equation  (40)  into  the  expres¬ 
sion  for  the  current  density, 


J  = 


B  J 


t)d3v 


(43) 


we  obtain  a  result  in  agreement  with  equations  (1)  and  (3) 


~  3  m  J  \ 


Vco86»t  4  ojsinMt  \  3  ®*o(v» 


2  2 
Cl)  4  v6 


•)v= 


8  v 


4  0(  Y  -2L)  +0(  y3)  . 
M 

(44) 


Thus  for  elastic  scattering  which  is  isotropic  in  the  center  of  mass  system 
we  have  proven  that  equation  (1)  rather  than  equation  (4)  is  the  correct 
expression  for  the  current  density;  however,  equation  (1)  is  valid  for  much 
more  general  conditions.  In  fact  it  is  valid  to  lowest  order  in  powers  of  y 
when  the  source  S(v,  t)  is  due  to  any  processes  provided  that  the  lowest 
order  anisotropy  in  S(v,  t)  is  of  order  or  higher  and  that  the  time  dependence 
of  S(v,  t)  is  of  order  Y  or  higher.  From  equation  (19)  we  see  that  the  latter 
condition  implies  that  the  purely  isotropic  part  of  S(v,  t)  is  equal  to  vfQ(v,  t) 
plus  terms  of  order  Y  or  higher.  Thus  the  lowest  order  contribution  to 
equation  (43)  for  the  conductivity  results  from  substitution  of  Vf0(v,  t)  for 
S(v,  t')  in  the  term  proportional  to  Y  in  equation  (17)  for  f(v,  t).  As  noted 
in  obtaining  equation  (40),  the  terms  proportional  to  (v/tvx  drop  out  in 
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the  integration  over  t1  and  with  the  aid  of  the  first  of  equations  (41)  one  ob- 

2 

tains  equation  (44)  with  corrections  of  order  y  or  higher  for  the  current 
density. 

In  order  to  demonstrate  clearly  the  velocity  dependence  of  fQ 
and  V  required  for  negative  conductivity,  we  integrate  equation  (44)  by 
parts,  obtaining 


j 


4g  e2  E  f  (  Vcosq  t  +  tsslntft)  y3f  ^  t- 
3  m  L  «2  +  v2 


(v,  t)dv 


r 

dv  L 


v3(  V  cos  cut  -f  josinjgt) 


w 


For  any  finite  u  the  first  term  must  vanish  in  order  that  the  energy  dencity 
of  the  electrons  be  finite.  We  take  the  scalar  product  of  j  with  the  electric 
field,  average  over  a  period  of  the  field,  and  obtain  the  average  power  given 
by  the  field  to  the  electrons 


3v2  V 


2(  V  + 


«2) 


+ 


v3(<o2  -  V2)  d V  \ 

2(  V2  ♦  a>2)2  dv  ' 


(45) 

In  the  vicinity  of  a  particular  speed  v^  we  approximate  the  collision  fre¬ 
quency  by  .V(vj)  *  avjn.  Then  it  is  apparent  that  for  ai>>V  (vj)  electrons 
of  speed  vj  lose  energy  to  the  field  (in  phase  part  of  the  conductivity  nega¬ 
tive)  if  n  <  -  3,  while  for  «<<  V(vj)  they  lose  energy  to  the  field  if  n>3. 
Thus  if  the  distribution  function  is  sharply  peaked  about  the  speed  Vj  for 
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which  V(v  j)  satisfies  either  of  these  conditions,  the  total  power  gained  by 
the  electrons  from  the  field  is  negative. 

It  is  possible  to  give  a  rough,  qualitative,  explanation  of  this 
negative  conductivity  phenomenon.  In  the  case  where  6>>>y  and  V  de¬ 
creases  with  velocity,  the  electron  oscillates  about  its  initial  velocity  many 
times  between  collisions  and  the  probability  that  it  will  undergo  a  collision 
during  the  portion  of  the  cycle  during  which  its  speed  is  below  its  initial 
speed  is  higher  than  the  probability  that  it  will  collide  during  the  portion  of 
cycle  when  its  speed  is  above  its  initial  speed.  It  will  thus  on  the  average 
lose  energy  to  the  field.  In  the  case  where  u«V  and  V  increases  with 
velocity  it  is  easiest  to  think  of  the  situation  as  approximating  the  D.  C. 

Case.  After  a  collision  an  electron  may  move  either  with  or  against  the 
field  (with  equal  probability  for  isotropic  scattering).  Since  the  collision 
frequency  increases  with  speed,  the  electrons  which  are  being  decelerated 
will  collide  later,  on  the  average,  than  those  that  are  being  accelerated. 

Thus  electrons  spend  more  time  giving  energy  up  to  the  field  than  gaining 
energy  from  it. 

Since  we  integrated  from  •«  to  t  in  obtaining  equations  (40)  and 
(44)  one  might  expect  our  equation  (44)  to  be  valid  only  when  fQ(v,  t)  is 
equal  to  the  steady  state  distribution  (the  Margenau  Distribution)  given  by 
equation  (37).  When  this  distribution  function  is  used  in  equation  (44) 
positive  values  result  for  the  conductivity,  regardless  of  the  velocity  dependence 
of  V.  Actually,  if  initially  the  electrons  have  some  strange,  non-isotropic, 
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distribution,  equations  (40)  and  (44)  become  applicable  after  a  time  equal 
to  only  a  few  multiples  of  the  collision  time  y"*.  Physically  the  time 
required  is  the  time  necessary  for  nearly  all  electrons  to  have  suffered 
one  or  more  collision,  mathematically  it  is  the  time  t  required  for  transient 
terms  proportional  to  exp  (-  V  t)  to  become  negligible.  This  time  is  really 
very  small  in  comparison  with  the  time  required  to  reach  a  steady  state. 

A  steady  state  is  reached  when  the  average  energy  gained  by  an  electron 
from  the  field  between  successive  collisions  is  equal  to  the  average  energy 
lost  by  an  electron  to  a  molecule  per  collision.  Both  these  energies  are 
small  compared  with  the  average  kinetic  energy  per  electron.  The  condi¬ 
tion  that  the  energy  gained  between  collisions  be  small  compared  with  the 
average  kinetic  energy  per  electron  must  bo  met  in  order  for  our  approach 
to  be  valid  (a  condition  met  by  most  physical  cases).  The  ratio  of  energy 
loss  per  collision  to  electron  kinetic  energy  is  ~m/M  ~10-^.  Thus  we 
expect  that  a  time  equal  to  thousands  of  collision  times  may  be  required 
before  a  steady  state  Is  reached.  As  an  example  suppose  that  Initially 
a  neutral  gas  in  a  fairly  weak  electric  field  is  bombarded  with  a  beam  of 
low  energy  electrons.  After  a  few  collision  times  the  electron  distribution 
is  changed  from  an  approximate  delta  function  of  velocity  to  an  approximate 
delta  function  of  energy  and  equations  (40),  (44)  and  (45)  apply.  If  the 
velocity  dependence  of  V(v)  happens  to  meet  the  conditions  discussed 
below  equation  (45),  energy  will  be  lost  by  the  electrons  to  the  field,  per¬ 
haps  for  a  detectable  length  of  time  if  the  density  of  the  gas  is  low  (  V  small). 


Z4 


© 


As  the  steady  state  is  approached  the  distribution  becomes  peaked  about 
the  origin  in  velocity  space  and  the  conductivity  becomes  positive  as 
seen  from  equation  (44). 

We  turn  now  to  a  consideration  of  the  contribution  of  higher  order 
terms  in  powers  of  Y  .  These  terms  might  be  expected  to  give  significant 
contributions  to  equations  (37),  (40)  and  (44)  in  the  case  of  strong  fields. 

In  this  connection  we  will  also  discuss  more  fully  the  criteria  for  the  vali¬ 
dity  of  the  expansion  in  powers  of  Y  UBed  in  obtaining  these  equations. 

In  principle  our  technique  can  be  used  to  obtain  the  part  of  f(v,  t) 

3 

which  is  proportional  to  Y  ;  however,  it  would  be  rather  tedious.  Instead 
we  show  by  a  more  careful  consideration  of  equation  (37)  that  such  a  calcu¬ 
lation  would  be  of  little  value  because  for  steady  state  conditions  the  con- 
tribution  of  the  next  higher  order  terms  in  powers  of  Y  to  either  fc(v,  t) 
or  j  are  of  order  m/M  or  less  compared  to  those  retained  in  equations 
(37)  and  (44). 

We  see  this  in  the  following  way:  if  the  second  term  in,  the  denomi¬ 
nator  of  the  exponential  in  equation  (37)  is  large  compared  with  the  first 
we  have 


•• 


MY 


,,  2 
6(  w  + 
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(46) 


or 
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(47) 
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2  2  2 
(  o>  +  V  )v 


where  Te  is  the  effective  electron  temperature  obtained  by  setting  kTe 
equal  to  the  denominator  of  the  exponent  of  equation  (37)  and  V  is  the  approxi¬ 
mate  average  magnitude  of  v.  On  the  other  hand,  if  kT  >  y^M(o)2  +  V2)”' 
2  -2  2  2  1 

we  see  that  y  (7)  (  u  4  V  )  is  even  smaller.  The  criterion  for 

the  validity  of  equation  (40)  is,  of  course,  that  the  ratio  of  successive  terms 
be  less  than  unity.  From  inspection  of  equations  (40)  and  (41)  this  criterion 
is  roughly 


Y 


w? 


<  1  ;  «>  V 


and 


<4B) 


- —  <  1  ,  »<  V 

W 

comparing  this  with  equation  (47)  we  see  that  in  the  steady  state  condition 
a  slightly  ionized  gas  adjusts  itself  automatically  in  such  a  way  that  the 
expansion  given  by  equation  (40)  is  valid  and  the  higher  order  corrections 
to  equations  (37)  and  (44)  are  of  order  m/M  or  less  relative  to  those  re¬ 
tained.  Of  course,  for  high  enough  field  intensities  kTe  becomes  so  large 
that  various  kinds  of  inelastic  and  ionizing  collisions  become  frequent. 
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Then  our  assumption  that  all  collisions  other  than  elastic  collisions  be¬ 
tween  electrons  and  molecules  are  insignificant  in  determining  f(v,  t) 
breaks  down.  Thus  the  only  circumstances  for  which  terms  of  order  y 
could  give  a  significant  contribution  to  equation  (44)  would  be  for  short 
times  after  an  electric  field  was  applied  to  a  Blightly  ionised  gas  in  which 
the  electrons  initially  had  energies  of  the  same  order  or  less  than  that 
gained  on  the  average  by  an  electron  from  the  field  between  successive 
collisions.  For  constant  collision  frequency  and  isotropic  distribution 
of  electrons  the  latter  energy  is  (  +  V2)"^  y2m/2.  Thus,  it  is  interest¬ 

ing  to  note  that  the  criteria  given  by  the  inequalities  (48)  for  the  validity 
of  our  expansion  in  powers  of  y  is  roughly  equivalent  to  demanding  that 
the  average  energy  gained  by  an  electron  from  the  field  between  collisions 
(assuming  V  Is  not  a  very  strong  function  of  v  )  be  small  compared  with 
the  average  kinetic  energy  per  electron. 

It  is  desirable  to  have  the  results  for  the  case  of  non-isotropic 
scattering  in  the  center  of  masB  system  as  well.  The  conventional  result 
for  this  case  is  that  equation  (37)  and  (44)  apply  with  V  replaced  with  V^, 
the  collision  frequency  for  momentum  transfer,  This  problem  is  consi¬ 
dered  in  the  appendix.  There  it  is  shown  what  conditions  must  be  met 
by  the  part  of  f(v,  t)  proportional  to  Legendre  Polynomials  P^  (vx/v)  with 
/  2  2  in  order  to  satisfy  the  requirements  for  the  validity  of  the  conventional 
result.  It  is  plausible  that  these  conditions  are  met,  at  least  for  the  steady 
state;  however,  in  contrast  to  the  case  of  isotropic  scattering,  we  are  unable 
to  show  rigorously  that  they  •’.  e  met. 
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APPENDIX 


It  is  our  purpose  here  to  generalize  to  the  case  of  non-isotropic 
scattering.  In  this  case  we  define  a  pseudo  source  function  Sm(v,  t)  by 
the  equation 


-  Vf(v,  t)  +  S(v,  t)  =  -  Vmf(v,  t)  +  Sm(v,  t) 


(49) 


where  is  the  collision  frequency  for  momentum  transfer 


Vm<v)  =  Nv 


(1  -  cos  ff)  a  (v,  0  )dfl  . 


(50) 


When  equation  (49)  is  used  to  express  the  right  hand  side  of  the  Boltzmann 
equation,  equation  (6),  we  obtain  the  generalized  versions  of  equations  (12), 
(17)  and  (18).  By  the  generalized  version  of  an  equation  we  will  mean  that 
the  only  change  from  the  original  is  that  V(v)  and  S(v,  t)  are  everywhere 
replaced  with  Vm(v)  and  Sm(v,  t),  The  symmetry  properties  of  Sm(v,  t) 
can  be  investigated  by  solving  equation  (49)  for  Sm  and  by  using  equations 
(25)  and  (27)  together  with  the  definition  of  V,  second  of  equations  (29). 

We  get 


Sm<I'  *> 


/n  [  fo(v,’t)  +  fl<v’-t> 


(V)4  a(v',0) dO 


-  Nv  [  f0(v,  t)  4  fjjv,  O  (v, $  )  cosfldfl 
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(51) 
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The  only  agency  producing  asymmetries  in  the  distribution  function  is  the 
electric  field.  Thus  the  purely  anisotropic  part  of  the  distribution  function 
can  be  expanded  in  Legendre  polynomials  as  follows: 


fj(v,  t)*  =  Cg[v,  t)  ¥  ff  (coster);  cos  *  »  v„/v  . 


crv,t)7 


(52) 


Now  we  substitute  this  expression  into  equation  (51)  and  expand  functions 
of  v' in  powers  of  m/M  about  v  with  the  use  of  equation  (26).  We  obtain 


S_(v,  t)  =  V  f 
m  m  o 


(v,  t)  +  Nv  2  C £  <v,  [7  (co8ii")  -  p p  (cos(fr)cos6^j  o (v,  fl)dO 

*  fl  (53) 

+  "  “  {  v3  VmfQ(v,  t)  +  Nv4  ^2  C/  <V)  (cos^'Ml- cos0)  ff(v,  0)dftj, 


We  now  make  use  of  the  well  known  relation 


Pa  (cos  \!>')  =  Pi  (cosy^)  P/ (co»  $ )  +  2  2  ^  ~  m> !  Pi.m(coBfA)  Pi1V'(cosfl)cos  m( ^  -  <k,) 

/  /  /  m=l  (/  +  m) !  /  2 

*  (54) 

where  (  £  -  ^2)  is  t*'e  an8^e  between  the  projections  of  v1  and  E  in  the  plane 
perpendicular  to  v.  Since  there  is  no  azimuthal  asymmetry  in  our  problem, 
the  part  of  equation  (54)  proportional  to  cosm(  £  -  £2)  gives  no  contribu¬ 
tion  when  equation  (54)  is  substituted  into  equation  (53).  Thus  we  obtain 
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Sm<I«  «  "  *mfo(v*  ‘>  +  f£z  Cf  (v*  *>  P/  <cos  *>  [  Vm<v>-  ^  (v)  ]  +  0(J|) 


in  which 


J(v)  =Nv^  £l  -  P^  (cosg)  j  o(v,  0  )d0 


In  contrast  to  this  result  we  recall  that  the  source  function  for  the  case 
of  isotropic  scattering  in  the  center  of  mass  system  had  no  zero  order  con¬ 
tribution  in  powers  of  m/M  arising  from  f  j(v,  t),  the  anisotropic  part  of  the 
distribution  function.  By  our  particular  choice  of  definition  for  the  pseudo 

source  function  S  {v,  t)  we  have  eliminated  the  zero  order  contribution 
m  ~ 

proportional  to  Pj(cos0)  =  costf».  However,  one  sees  that  if  Sm(v,  t) 
had  been  defined  by  equations  (49)  and  (50)  with  any  other  Legendre  poly¬ 
nomial  P£  (cos0)  replacing  cos  $  in  equation  (50),  then  the  zero  order 
anisotropy  proportional  to  Pj  (cos' l>)  would  have  been  eliminated  from 
S_(v,  t).  Nevertheless  the  choice  made  for  the  definition  of  S^(v,  t)  is 
most  convenient. 

If  we  make  the  assumption  (usually  made)  that  the  part  of  f(v,  t) 
proportional  to  Legendre  polynomials  (cos^)  with  /  <S  2  is  negligible, 
S^^v,  t)  is  sufficiently  isotropic  for  us  to  proceed  as  in  Section  11,  We 
then  get  the  conventional  result  that  the  steady  state  value  of  f0(v,  t)  is 
given  by  the  generalized  version  of  equation  (37)  and  the  current  density 
by  the  generalized  versions  of  equations  (1)  or  (44).  Actually  all  that  is 
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required  for  these  results  to  apply  is  that  the  C ^  (v,  t)  with  /  2  2  be  of 
first  order  or  higher  in  powers  of  y,  When  the  generalized  version  of 
equation  (17)  is  integrated  over  angle  to  get  the  generalized  versions  of 
equations  (19)  and  (20),  the  (v,  t)  with  f  £  2  appear  only  in  terms  pro¬ 
portional  to  y2  or  m/M  as  a  consequence  of  the  orthoganality  properties 
of  Legendre  polynomials.  Thus,  if  all  C p  (v,  t)  with  !  2  2  are  of  order 
y  or  higher,  they  contribute  terms  to  equations  (19)  and  (20)  which  are 
of  order  y *  and  ym/M  or  higher  and  can  be  neglected.  We  will  show 
now  that  the  requirement  that  the  C ^  (v,  t)  be  of  order  Y  or  higher  for 
f  2  2  reduces  to  the  requirement  that!  1  -  Vm/iy  )  times  the  approximate  average 
change  in  C ^  (v,  t)  between  collisions  be  of  order  V  or  higher.  An  equa¬ 
tion  for  C2(v.  t)  is  obtained  by  substituting  the  value  for  the  pseudo  source 
function  given  by  equation  (55)  into  the  generalized  version  of  equation  (17), 
Multiplying  the  equation  through  by  P^tcos^),  and  integrating  over  solid 
angle,  one  gets 


C2(v'  l)  *  ( Vm  ‘  V2)  /  e*  Vm<t  "  l,)  C2(v’  t,)dt'  +  °<  V  2> 


which  can  be  rearranged  as 


c2(v.  t)  =  ( 1  -  -t~)  < 4  C2(v>  >  +  °<  y 2) 
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( 


<  4C2(v,  t)  >  is  approximately  the  average  change  in  C2(v,  t)  betv/een 
successive  collisions  given  by 


V 

<AC2(v,  t)>  =  /  Vmdt'  [c2(v,  t)-C2(v,  t')]  . 


Proceeding  in  an  analogous  fashion  we  can  get  similar  expressions  for 
the  Cj  (v,  t)  with  H>  2.  For  j!  >  Z  the  analogues  to  equations  (57)  and 
(58)  have  a  term  going  as  y  or  higher  in  place  of  the  y  term.  Thus, 
in  all  cases  the  requirement  that  the  C£  (v,  t)  with  ft>  Z  be  of  order  y 
or  higher,  which  is  necessary  for  the  validity  of  the  generalized  ver¬ 
sions  of  equations  (37),  (1)  and  (44),  is  seen  to  reduce  to  the  requirement 
that  ( 1  -  Vm/  Vj^)  <  A  Cj  ( v,  t)  >  be  of  order  Y  or  higher  for  all  jP  t  Z, 

In  the  fairly  low  energy  range  where  one  is  Justified  in  ignor¬ 
ing  inelastic  processes,  the  DeBroigle  wave  of  the  electron  can  be  expected 
to  be  long  compared  to  the  range  of  molecular  forceB,  in  which  case  the 
scattering  will  be  mostly  s-wave  scattering.  Then  S  Vm  and 

the  first  term  on  the  r.h.  s.  of  equation  (58)  will  be  small,  (This  would 
not  hold  for  molecules  with  appreciable  static  dipole  moments.)  However, 
separately  we  expect  <ACg  (v,  t)>  to  be  small,  particularly  for  the 
steady  state.  When  a  steady  state  is  obtained,  we  could  separate  (v,  t) 
into  a  constant  part,  which  would  contribute  nothing  to  <ACg[v,  t)  >  , 
and  an  oscillating  part.  Since  the  only  agency  tending  to  produce  oscil¬ 
latory  motion  is  the  electric  field,  we  would  expect  that  the  oscillations 
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would  be  at  one  or  more  harmonics  of  the  field  frequency  and  that  this 
part  ol  C  a  (v,  t)  and  hence  <AC^(v,  t)  >  ,  would  go  as  some  power  of 
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